Universality in the jamming limit for elongated hard particles in 
one dimension 
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Abstract. - We study thermodynamics properties of a one dimensional gas of hard elongated 
particles. The particle centers are restricted to a line, while they can rotate in two-dimensional 
space. Correlations between orientations of the objects are studied (by transfer matrix method) 
as a function of density and aspect ratio. The behavior in the extreme high-density (jamming) 
limit is described by a few universality classes depending on the object's shape. In particular, 
there is a diverging correlation length when the contact point of adjacent objects is far from the 
line along which their centers move, as for needles and rectangles. 



One-dimensional (ID) collections of classical particles 
are usually simple to study, yet can offer valuable insights 
into collective properties of higher dimensional systems. In 
particular, thermodynamic properties of a gas of particles 
confined to a line, and interacting with a short-range po- 
tential can be solved exactly [1]. For example, the solvable 
ID system of "hard spheres," sometimes referred to as the 
Tonks gas [2] , served as an initial step in the study of two- 
and three-dimensional systems of hard disks/spheres, and 
more recently was invoked in connection with proteins on 
DNA [3]. 

Usage of hard particles, whose interaction energy is ei- 
ther or oo, is yet another useful simplification that high- 
lights the geometric/entropic features of interacting sys- 
tems [4]. The absence of an energy scale renders the be- 
havior independent of temperature, emphasizing the varia- 
tions with shape and density. Indeed, numerical studies of 
hard potentials date back to the origins of the Metropolis 
Monte Carlo method [5]. More recently, dynamic proper- 
ties of hard particles have been studied in connection with 
granular flows [6] and jamming transitions [7]. Critical be- 
havior and diverging correlation lengths at the so-called 
J-point (for jamming) have been studied from several per- 
spectives [8-11]. 

Here, we demonstrate that diverging correlation lengths 
may also arise in the high density (jamming) limit of elon- 
gated hard particles in one-dimension (such as depicted in 
Fig. 1). Non-spherically symmetric hard particles [12-24] 



provide important perspectives into the phase behavior of 
liquid crystals. Indeed, following an approach outlined in 
Ref. [25], Lebowitz et al. [26] studied this very system of 
elongated objects moving on a line and rotating in the 
plane. They showed that in the limit of high density, the 
ID pressure acquires a simple form, dependent on the cur- 
vatures of the object at their point of contact. We build 
upon these results, focusing on the orientational fluctua- 
tions and their correlations. We show that, depending on 
the shapes of the objects, there are a few generic univer- 
sality classes in the jamming limit, some characterized by 
diverging correlation lengths. 

As indicated in Fig. 1, the elongated particles consid- 
ered in our work are characterized by the positions {xi\ 
of their centers of mass, and their orientations {<j)i\ (mea- 
sured with respect to the normal to the line) . We consider 
symmetric objects that are invariant under a change of 
(j)i by TT, and thus restrict — 7r/2 < (pi < tt/2. (There are 
obvious generalizations to three dimensional objects, with 
more angular degrees of freedom, which shall not be pur- 
sued here.) The elongated objects are typically described 
by long and short semi-axis lengths a and b {a > b). The 
needle can be viewed as either an ellipse or a rectangle 
with 6 = 0. The system consists of N identical particles, 
and without loss of generality we shall assume that their 
order along the line remains unchanged, i.e. Xi < x^+i. 

The particles are not allowed to intersect but do 
not interact otherwise. The minimal distance di^i i = 
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Fig. 1: (Color online) Some elongated objects considered 
in this work: (a) needles, (b) rectangles, (c) ellipses, and 
(d) metronomes. Centers of these two-dimensional objects are 
placed on a straight line. The degrees of freedom for each par- 
ticle are the position of its center Xi (constrained to ID), and 
the angle <j)i between its major axis and the normal to the line, 
as depicted in (a). The longer and shorter semi-axes of the 
particles are denoted a and b, respectively. All particles in a 
particular ensemble are identical. Needles can be viewed as the 
b ^ limit of either ellipses or rectangles. 



d{(pi-i, between two adjacent particles depends on 
their orientations (f>i-i and and the size and shape 
of the objects. The hard particle constraint is equivalent 
to the interaction 



0, if Xi - Xi^i > di-i^i, 
oo, otherwise. 



(1) 



The standard method for dealing with potentials that 
depend only on the nearest neighbor separations is to con- 
sider an ensemble of fixed pressure p, and evaluate the 
Gibbs partition function [27, 28] at temperature T: 
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where /? = l/ksT, and fc^ is the Boltzmann constant. In 
the second line of Eq. (2), we have changed variables to 
the separations setting the boundary con- 

ditions xq = 4>o = for convenience. Since the potential 



V is zero or infinity, the integration over each Ui is easily 
performed, leading to the final result, which only depends 
on the combination Pp. 

As pressure is increased the particles are jammed to- 
gether, and forced to align to minimize their separations. 
The thermodynamic properties of the system in this limit 
can be obtained by evaluating the Gibbs partition function 
in Eq. (2) by a saddle point expansion around {(pi = 0}. 
The expansion of d{4>, (f>') depends on the shape of the ob- 
jects, and for the cases depicted in Fig. 1 we have: 



• For rectangles 



(3) 



The first two terms reflect the change in the cross- 
section of each rectangle as it tilts. The flnal term 
reflects the advantage of aligning neighboring rectan- 
gles. 

• Needles are obtained as a limit of the above expression 
when 6 = 0. The full expression, including higher 
order terms is [29] 



d(0,0') 
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• The minimal distance between two ellipses can be ob- 
tained by using the expressions in Ref. [30] . For small 
tilts from vertical, we find 
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Note that while ellipses become needles as 6 0, 
Eq. (4) corresponds to a different branch of solutions, 
and is not obtained as a limit of the above expres- 
sion. Equation (5) thus cannot be used to examine 
the crossover to needles. 

• Metronomes are a singular limit in which d — 2b for 
small tilts (and given by Eq. (4) for larger tilts) . The 
high density configuration is not unique and resembles 
that of hard circles. 

From these results, we observe that the general expan- 
sion for symmetric objects has the form 



d(0, 4>') =do + di{<j>^ + 0'2) -I- d2 10 - 0T + 



(6) 



Here, do indicates the separation at jamming, while di 
indicates how the footprint of the particle is increased 
as it tilts. For needles do = di — 0, while c?i = for 
metronomes. The correlations between orientations are 
due to the third term in the expansion, and we note that 
generically 7=1 when the intersection point is away from 
the line (as in needles and rectangles), while 7 = 2 when it 
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is close to zero for small tilts (as for ellipsoids). Lebowitz 
et al. [26] allow for the possibility of a general 7, although 
this will only be the case for singular shapes and is not 
generic. In any case, for 7 < 2 we can ignore the sub- 
dominant term proportional to di in a saddle-point evalu- 
ation of the final integrals in Eq. (2). We can then change 
variables to it) = f3pd2 \(f> — (fi'f, and simple dimensional 
analysis leads to Z cx e-^^P*(/3p)-^(i+''"'), and 
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We can define the free space available for fluctuating par- 
ticles as « (1 -l-7~^)//3p, and in the remainder we shall 
use this expression to switch between (jp and £/. (Note 
that the saddle point estimate is only asymptotically exact 
at high pressures.) 

To study orientational fluctuations and correlations, we 
coarse-grain the variables {(jii} to a continuous field (/)(n). 
From symmetry grounds, we expect that configurations of 
the coarse-grained field are distributed according to the 
Boltzmann weight of a Hamiltonian 
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For the case of ellipses, with 7 = 2, Eqs. (2,6) imme- 
diately suggest J ~ 4/3pdi and (replacing differences by 
derivatives) K « 2/3pd2. In the limit of a ^ 6, Eq. (5) 
leads to J « 2/3p6 and K « /3paV(26). 

It is less clear how a discrete sum of distances as in 
Eq. (4) for needles goes over to the continuum form of 
Eq. (8). However, if we integrate (in the saddle point 
sense) the weight in Eq. (2) over every other (say odd) (f)i , 
the effective weight over the remaining (even) 4>iS takes 
the form 



i odd 



[{(3paf{(j>i+i - (p^ 



(9) 



In the continuum limit this then leads to J « 1/2 and 
K « 2{l3paY . The appearance of a finite J « 1/2 is 
worth noting: while at the lowest order there is no cost to 
tilting a needle, the nonlinearities implicit in the denomi- 
nator of Eq. (4) create a preference for = after a single 
integration. While for needles this is the leading term, in 
general we expect such contributions from non-linearities 
as corrections to the larger terms proportional to powers of 
fip. To recover the above result for needles as (6/a) 0, 
we may employ the interpolation formula J ~ 2j3ph +1/2 
for both ellipses and rectangles (from Eq. (3)). While we 
shall use this interpolation form in the remainder of the 
paper, we do not expect it to hold exactly (discreteness 
and non-linearities are likely to lead to a more compli- 
cated formula). As in needles, we get K w 2{fipaf for 
rectangles, while in the case of metronomes, we expect a 
finite value of K to result from non-linearities. 




0.001 



/ 10 



Fig. 2: (Color online) Variance of angle as a function of dimen- 
sionless pressure / = Ppa for rectangles. The rightmost (red) 
solid line is for needles (a = 0), while the remaining solid lines 
correspond to aspect ratios a — 0.01,0.02,0.05,0.1,0.2,0.5 
(right-to-left). The lower dotted line demonstrates the fast 
decay for a = 0.7, while the almost horizontal dotted line cor- 
responds to a square particle q = 1. Short straight segments 
of slopes -1 and -1.5 are included for visual comparison. 



From the Gaussian weight associated with Eq. (8), it 
is easily obtained that the fluctuations in angle for each 
object are characterized by a variance cr^ = (0^) = 
(4JA')~^/^. Using the asymptotic forms for J and K ob- 
tained above, this leads to 



2/\/l+4a/' 
^/(l+4a/) ' 



for rectangles 
for ellipses 



(10) 



where we have introduced the dimensionless parameters - 
reduced pressure / — Ppa and aspect ratio a — b/a. We 
numerically checked these results by utilizing the trans- 
fer matrix [25,28,31] M{(t),(j)') = g-'^P'^^"*'"*') to evaluate 
the integrals in Eq. (2). The angular variables were dis- 
cretized to obtain finite matrices, which were then multi- 
plied together as detailed in Ref. [29], where the transfer 
matrix method was employed for computation of elastic 
moduli. Naturally in the numerical computations we used 
the exact form of the separation d{<j), (j>') for rectangles and 
ellipses, and not the asymptotic forms given above for the 
saddle point evaluations. 

Figure 2 depicts the variance in angular fluctuations, as 
a function of reduced pressure /, for rectangles with sev- 
eral aspect ratios. For / < 1, all graphs have a constant 
value corresponding to a uniform distribution of orienta- 
tions, since in this dilute gas regime the orientations are 
unconstrained. For larger /, the curves for different val- 
ues of a initially follow the behavior of a needle (a = 0), 
with (T^ cx 1//, but then depart, and eventually scale as 
~ l//'^^^i which characterizes the variance of rectan- 
gles. The crossover occurs at larger / for smaller aspect 
ratios, as predicted by Eq. (10). The case of a = 1 is 
exceptional, because the rectangle becomes a square, with 
two preferred orientations instead of one. 
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Fig. 3: (Color online) Variance of angle as a function of dimen- 
sionless pressure / for ellipses. The rightmost (red) solid line is 
for needles (q = 0), while remaining solid lines correspond to 
aspect ratios a = 0.01, 0.02, 0.05, 0.1, 0.2 (right-to-left). Dotted 
lines (left-to-right) demonstrate the function for a — 0.5, 0.7, 
while the horizontal dotted line correspond to a circular parti- 
cle a = 1. The short straight segment of slope -1 is included 
for visual comparison. 



The pressure-dependence of the variance of the angle for 
ellipses with different aspect ratios is depicted in Fig. 3. 
Again all curves coincide for / < 1, with a value cor- 
responding to a uniform distribution in angles. Upon 
increasing /, the variance approaches the scaling form 
o"^ ^ 1// predicted by Eq. (10). Note that the overaU 
trend of the graphs moving right-to-left with increasing a 
is reversed for a = 0.5, and eventually for a = 1 the ellipse 
becomes a circle with all orientations equally probable, in- 
dependent of /. 

We also numerically examined correlations of the ne- 
matic order parameter s = cos 20. The two-point corre- 
lation function c{n) = {s{i)s{i + n)) — (s)^ decays from 
(cos^(2</))) for n = to as n — > oo. From the Gaussian 
weight associated with Eq. (8), it can be verified that 



c(n) = 2 exp 



sinh 



(11) 



Specializing further to the limit / ^ 1, required for the 
validity of the Gaussian approximation, we find c(n) cx 



-n/5 



, with ^ = K/{4:J), and in particular 
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for rectangles 
for ellipses 



(12) 



Note that the correlation length diverges as / ^ oo for 
rectangles, but asymptotes to a constant for ellipses. 

Correlation functions were obtained numerically, by ap- 
propriate transfer matrix manipulations, for a variety of 
rectangles and ellipses with different aspect ratios, and 
confirm the above predictions. In particular, we verified 
Eqs. (12) for a moderate value of the aspect ratio [a — 0.05 
for rectangles, and a = 0.1 for ellipses) chosen such that 
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Fig. 4: (Color online) Numerical verification of Eqs. (12) for 
the asymptotic relation between the correlation length ^ and 
the dimensionless pressure /. The top panel is for rectangles 
with aspect ratio a = 0.05, and the bottom panel is for ellipses 
with aspect ratio a — 0.1. 



the crossovers in / are clearly seen. In the plots of Fig. 4, 
^ is rescaled in such a way that Eq. (12) is verified if the 
curve is a straight line. Indeed for large values of /, where 
the Gaussian approximation is valid: the dependence is a 
straight line. The slopes and the intercept points of the 
lines arc of the order predicted by Eq. (12) but do not 
agree quantitatively, since the approximate nature of our 
derivation does not ensure accuracy of the coefficients in 
Eq. (12). 

The most notable feature of our results is the possible 
divergence of the dimensionless correlation length ^. This 
diverging correlation length cannot be associated with the 
onset of long-range order (as in say the ID Ising model) , as 

there is a finite nematic order (cos(20)) « exp 

at any density. As such, the onset of jamming [7] is a bet- 
ter descriptor: As pressure/density increase, and the free 
space £f vanishes, a larger collection of particles must be 
moved in concert. This number (correlation 'length') di- 
verges for rectangles but saturates for ellipses (see Eq. (12) 
for / — > oo). The differences between the two cases is 
that for rectangles the point of contact is remote from 
the center, almost at a distance a. Thus, given a typical 
spacing of near-neighbor angular shifts are small and 
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vanish as ~ ^//a. By contrast, for ellipses the point 
of contact goes to zero continuously, roughly as ^ Rl j , 
where R is the radius of curvature at the point of con- 
tact. This results in larger near-neighbor angular shifts 
of Sep ^ \JljjR. At the global scale, fluctuations in angle 
are limited (due to the broken nematic symmetry imposed 
by the line of movement) and given by Eq. (10), accord- 
ing to which tJ is proportional to (^'j'^ , and for 
needles, rectangles, and ellipses respectively. As If ^ 0, 
the local (near-neighbor) fluctuations are asymptotically 
smaller than global ones for rectangles (and needles) but 
not for ellipses. The fluctuations must thus build over a 
large correlation number in the former cases. 
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